%不同情况，需要更改调用的函数
%龙格库塔函数  解微分方程
% [0:0.1:20] t的范围，0.1是步长
[t,w]=ode23(@ifun1,[0:0.1:20],[0,0.01,0,6,0,0.01],[ ]);
plot3(w(:,1),w(:,3),w(:,5));%绘制轨迹
grid on
title('b\neq0,e=0')   %图片标题